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Abstract. The phase diagram of cold dense quark matter is studied using the’t Hooft extended Nambu-Jona-Lasinio Model 
applied to the light quark sector with a finite current mass for the strange quark (up and down are considered in the chiral 
limit). By relaxing the traditional uniformity assumptions and considering a modulated light quark condensate background 
we investigate the possible existence of non-uniform phases in this region of the phase diagram. The effects of changes in the 
coupling strengths of the model are studied and it is shown that the inclusion of flavour mixing combined with the finite current 
mass of the strange quark catalyses the appearance of the non-uniform phases, extending the domain for their existence. 
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INTRODUCTION 

Low energy models of the Nambu-Jona-Lasinio (NJL) [1, 2, 3] type have long been used as a tool in the theoretical 
study of strongly interacting matter. The success of this approach can be tied to the fact that they share with QCD 
its global symmetries and incorporate a mechanism for the dynamical breaking of chiral symmetry. Explicit breaking 
of f//i(l) can be incorporated through the inclusion of the OZI-violating ’t Hooft flavour determinant [4, 5, 6 , 7, 8 ]. 
These simple models are of particular usefulness in the low temperature/high chemical potential regime due to the sign 
problem which affects the more first principle approach of lattice-QCD. 

As the pionic interaction with quarks (or nucleons) is attractive when the mean pion field carries a gradient it has 
long ago been proposed (for recent reviews see, for instance, [9, 10]) that the non-trivial balance between this effect 
and that of the kinetic terms could allow for the existence of non-uniform phases in cold dense strongly-interacting 
matter. 

Here we will outline some of the results reported in [11] (to which we refer for a more in-depth discussion), 
pertaining to the case study of dense light-quark matter (we consider the u, d and s quarks) at zero temperature 
using the NJL model with a’t Hooft determinant term. A non-vanishing current mass for the s quark is taken into 
account whereas the chiral limit is considered for the u and d quarks thus validating the use of a simple chiral wave 
ansatz for the light quark chiral condensates. 


THE MODEL 

In the chiral limit the Euler-Lagrange equations of the model admit as an analytical solution the chiral wave ansatz 
proposed in [ 12 ] and given by the following quark orbitals and energy levels: 

(Wim) = ycos(^-r), (W*75T3r;) = ysin(^-r), 

-\- p^ ± \J{p ■ qf' (1) 

Here M and p denote the dynamical mass and the momentum of the quark, q is the wave vector of the condensate 
modulation and T 3 is the Pauli matrix acting in isospin space. Note that the branch has a lower energy than £+ 



and thus this branch is preferably occupied. For the s quark no modulation is considered (uniform strange condensate 
background). 

In the mean field approximation the application of techniques of Ref. [13] yields the thermodynamic potential of 
the model as: 
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where /i, (i = u,d,s) are twice the quark condensates. The coupling strengths G and K refer respectively to the NJL and 
the’t Hooft determinant terms. The bosonization of the model through the introduction of auxiliary variables enables 
the evaluation of the fermionic part which is rendered quadratic. The stationary phase contribution coming from the 
integration over the auxiliary fields is given by while the contribution coming from the fermionic path integral is 
given by the integrals. The notation refers to the subtraction of the quantities evaluated at the dynamical masses 
by their value evaluated at M = 0 [14]. Using a regularization kernel corresponding to two Pauli-Villars subtractions 
in the integrand [15], previously used for instance in [16, 17], namely p (sA^) = 1 — (1 +sA^)exp(—sA^), the vacuum 
and medium contributions, and can be written explicitly as: 


T rVac I jMed 

7-1 =J_\ +7_i , 



where = y + A^. The z-axis was chosen to coincide with the direction of modulation wave vector. The 

notation refers to the subtraction of the same quantity evaluated for a uniform gas of massless quarks (M — 0 and 
q — 0). The superscript and subscript ± in the definition of iF refer respectively to the energy branch and the sign 
in front of the chemical potential in the exponent. The C(7’,p) term is added for thermodynamic consistency [14] 
ensuring that the correct asymptotic behaviour is obtained in the T —> oo limit. 

The minimization of the thermodynamical potential with respect to M and q has to be done self-consistently with 
the resolution of the stationary phase equations: 

{ Tiiu Mu = Ghu -\~ 

nid - Md = Ghd + ^Khs , (4) 

TTls Ms = Ghs -\~ J^^u^d 

where m, stands for the current masses of the quarks. 

RESULTS 

In our study the current masses were set to niu = nid = 0, nis = 186 MeV. The remaining three parameters of the 
model (G, K and A) can be reduced to two, the ’t Hooft coupling strength, K, and the dimensionless curvature, 
T = NcGA^/{2n^), by fitting A to reproduce a reasonable value for the vacuum dynamical mass of the light quarks 
(Ml = 300 MeV). 

With no OZI-violating term, the light and strange sectors are decoupled and in the chiral limit T = 1 is the critical 
value above which dynamical chiral symmetry breaking occurs: as a crossover for 1 < T < 1.23 and as a first order 



transition for higher values. The transitions are independent for the light and strange sector, with the latter shifted to 
higher jj. due to the current s-quark mass (both occur at values of jj. close to the vacuum dynamical mass 

For high enough chemical potential an energetically favourable modulated solution always emerges. This solution 
goes asymptotically to lim^_i.oc {h,q} = {0,2/r} and becomes degenerate with the trivial one. 

Above 1.23 < T < 1.53, a window for energetically favourable finite-q' solutions appears in the vicinity of the first- 
order transition associated with the condensate of the light quarks and the usual transition to a vanishing condensate 
is superseded by a transition to a finite but modulated condensate: {/i/,0} —>^ The window of existence of this 

phase ends with the light condensate going continuously to zero (note that if the condensate vanishes its modulation is 
ill defined). The extent of this windows is increased for higher values of T and for T > 1.53 it merges with the one at 
higher chemical potentials. 
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FIGURE 1. From left to right: the chemical potential dependence of the zero temperature solutions for the light condensate, hi, 
the strange condensate, hs, and for the modulation wave vector q. From top to bottom: each row corresponds to a different choice 
of the’t Hooft coupling strength ([k"] = GeV^). Thicker lines denote the inhomogeneous solutions. In the grey regions the value of 
q is undetermined, since hi = 0. 

The inclusion of the’t Hooft determinant term couples the three flavour sectors leading to several different scenarios, 
at a fixed curvature, depending on the the value of its coupling strength. Some examples of this can be seen in 
Fig. 1 for the T = 1.4 case. For —K> 290 GeV^^ a new window of energetically favourable modulated solutions 
appears in the vicinity of /r ^ (see the shark-fin shape for hi in the first row of Fig. 1). There are therefore 





































three intervals of chemical potential where the global minimum of the thermodynamical potential corresponds to a 
modulated solution. The boundaries of these intervals corresopnd to two first order transitions and three-crossovers, 
with the latter corresponding to the (dis)appearance of a non-vanishing hi. The two first order transitions exclude the 
occurrence of the q — Q transitions as they occur at slightly lower/higher chemical potential for the one taking place 
near A zoom of the behaviour of the chiral condensates near the transitions can be seen in Fig. 2 for the 

— K = 500 GeV^^ case. 

By increasing the coupling of the’t Hooft interaction even further we can expand the extent of these chemical 
potential windows with finite-^' solutions and for —K> 935 GeV^^ the first two merge. As a consequence two 
crossover transitions disappear (see second row in Fig. 1). For —K> 1660 GeV^^ the last transition is no longer 
to homogeneous solution as there is an overlap of two distinct modulated solutions and the first-order transition occurs 
between two phases with finite q (third row of Fig. 1 and Fig. 3). 
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FIGURE 2. Zoom for the chemical potential dependence of the chiral condensate solutions, with thicker lines referring to the 
finite-^ solutions, in the vicinity of the transitions (marked by the vertical dotted lines), case of k: = —500 GeV^^. 

The main results can be summarized in the form of critical chemical potentials associated with the (dis)appearance 
of the inhomogeneous phases shown in Fig. 4. 



(from left to right), showing the merging of the solution branches for strong enough flavor mixing. 

We conclude that the flavour mixing effect of the ’t Hooft determinant combined with the inclusion of a finite current 
mass for the strange quark acts as a catalyst for the appearance of globally stable inhomogeneous solutions in cold 
quark matter by coupling the behaviour of the light and strange sectors. Further extensions of this model to include, 
for instance, the eight quark interactions term [18, 19] are straightforward and can provide an even richer picture for 
the phase diagram of cold strongly interacting matter. 
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FIGURE 4. In 4a we present the critical chemical potentials as a function of T in the NIL model and in 4b as functions of K 
([«■] = GeV^) in the NJLH model for the T = 1.4 case. The upper dotted line corresponds to the cut-off A. The chemical potential 
of the first order transitions are marked with the full black lines. Dashed lines indicate the borders of the region where the finite- 
q solutions exist. We distinguish between three types of critical chemical potentials and an example for this distinction in the 
K = —500 GeV^^ case appears in Fig. 4c (type I with finite q and hi, type II with finite hi and vanishing q, and type III with finite 
q but vanishing hi). 
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